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High-resolution spectroscopy of the low-mass X-ray binary 
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Abstract. We present initial results from observations of the low-mass X-ray binary EXO 0748—67 with the 
Reflection Grating Spectrometer on board the XMM-Newton Observatory. The spectra exhibit discrete structure 
due to absorption and emission from ionized neon, oxygen, and nitrogen. We use the quantitative constraints 
imposed by the spectral features to develop an empirical model of the circumsource material. This consists of a 
thickened accretion disk with emission and absorption in the plasma orbiting high above the binary plane. This 
model presents challenges to current theories of accretion in X-ray binary systems. 
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1. Introduction 

X-ray spectroscopic observations should provide a useful 
probe of the accretion processes in low-mass X-ray bina- 
ries (LMXB's). The continuum emission from the bright 
central source is reprocessed as the photons are absorbed 
and re-emitted in the surrounding material. The details of 
the resulting X-ray spectra are extremely sensitive to the 
physical conditions in the plasma and therefore provide an 
excellent way to constrain the circumstellar environment 
and the mass flow that fuels the X-ray emission. Previous 
observatories have not had the spectral resolution neces- 
sary to resolve the discrete structure in the X-ray spec- 
trum. The Reflection Grating Spectrometer (RGS) on 
the recently launched XMM-Newton Observatory provides 
both a high resolving power and a large effective area mak- 
ing it uniquely sensitive to the diagnostic spectral features, 
particularly in the soft X-ray band that contains transi- 
tions from highly ionized charge states of the most abun- 
dant elements. 

In this letter we present the results of RGS obser- 
vations of EXO 0748-67, a highly variable LMX B tha t 
was first discovered with EXOSAT (Parmar et al. 19861 ). 
Analysis of the EXOSAT light curves, which show deep 
eclipses with a 3.82 hour orbital period and a complex 
dipping structure, led to a derived inclination angle of 
75° - 82° (Parmar et al. |1986D. The detection of type I 
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X-ray bursts confirmed th at the compact object is a neu- 
tron star (Gottwald et al. 1986 ). Observations with both 
ASCA (Thomas et al. |l997| ) and ROSAT (Schulz [l999| ) 
revealed structure in the soft X-ray spectrum, but the 
limited spectral resolution of these instruments made it 
impossible to distinguish whether this was due to absorp- 
tion or emission features. In the RGS observations that we 
present below, we find a spectrum that is rich in discrete 
structure including absorption and emission from ionized 
oxygen, nitrogen, and neon. We use the available spectral 
diagnostics to construct a quantitative empirical model of 
the circumsource flow. 

2. Data Reduction 

The RGS covers the wavelength range of 5 to 35 A with 
a resolution of 0.05 A (roughly constant across the band) 
and a peak effective area of ^ 140 cm^ at 15 A. A com- 
plete description of the instrument can be found in den 
Herder et al. ( |2001| ). EXO 0748-67 was observed repeat- 
edly during both the comissioning and calibration phases 
of the mission. We have chosen the two longest exposures 
for this analysis; the first was 49.3 ks on 2000 March 
28, and the second was 44.8 ks on 2000 April 21. The 
raw data were processed with the development version of 
the XMM-Newton Science Analysis Software (SAS). The 
spectra were extracted by first applying a 30" wide spa- 
tial filter to the CCD image. The surviving events were 
then plotted in dispersion channel vs. CCD pulse-height 
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space and additional filters were applied to collect events 
from the different spectral orders. Background spectra 
were extracted by applying the same spectral order fil- 
ters to events from a spatial region that was offset from 
the source location on the CCD image. 

We assigned nominal wavelengths to each dispersion 
channel based on the geometry of the instrument and 
the pointing angles of the spacecraft. The pointings were 
found to be offset by 46" for the first observation and 60" 
for the second using images from the European Photon 
Imaging Camera (EPIC). We expect the wavelengths to 
be accurate to ^ 0.010 A. The development version of the 
SAS did not yet have the capability to correct the events 
for aspect variations. However, after excluding the initial 
part of each observation (8.6 ks of the first observation 
and 1.7 ks of the second), we find that the pointings are 
stable to less than ^ 2". Since this is much smaller than 
the resolution of the telescope, further aspect correction 
is unnecessary. 

We have determined the effective area for each expo- 
sure by applying the same extraction regions that were 
used with the data to the full response matrix, which in- 
cludes all information on the efficiency of the instruments. 
Based on our ground calibration we expect the uncertainty 
to be less than ^ 10% for wavelengths longer than 9 A and 
at most ^ 20% at the shortest wavelengths. We have cal- 
culated the fiux for each observation using these effective 
area curves and the standard SAS exposure maps, which 
give exposure times for each wavelength bin. We combined 
the first order spectra from both instruments for the two 
observations in order to maximize the statistical quality 
of our measurement. 



3. Observational Features 

3.1. Light Curve 

The observed EXO 0748—67 count rate varies by up to a 
factor of ten on timescales as short as a few hundred sec- 
onds. This is illustrated in Fig. |l|, which shows a portion of 
the RGS light curve. There are periods of stable low-level 
emission where the count rate is less than ^ 0.5cts~^, 
periods where the count rate varies rapidly between ~ 0.5 
and 3.5cts~^, and several type I X-ray bursts with peak 
intensities from 2.5 to 5.2 ct a~^. We see no correlation be- 
tween the variability and the orbital phase; similar to the 

|199|), 



soft X-ray light curve from ASCA (Church et al 
the RGS light curve does not show the quiescent level and 
dipping structure that is characteristic of the light curves 
from higher energy observations. Most importantly, we see 
no eclipses during the times predicted for the hard X-ray 
eclipses. 

3.2. Discrete Spectral Structure 

As shown in Fig. || the EXO 0748-67 spect ra show signif- 
icant discrete structure both in absorption and emission. 
We see bright emission lines from O viii Lya and the O vii 




Fig. 1. Partial first order light curve in 100 s bins. The 
background level is indicated by the thick line. A repre- 
sentative error bar (for statistical error only) is shown at 
^ 4hr. The horizontal bar marks the location of the hard 
X-ray eclipse. The time units have an arbitrary origin. 



He-like complex. The Ne x Lya, Ne IX He-like complex and 
the Nvii Lya emission lines are weaker, but clearly visi- 
ble, particularly during the periods of low emission when 
the equivalent widths of the lines are highest. We see the 
photoelectric absorption edges of both O viii and O vii, 
particularly during the periods of rapid variation when 
the contiuum intensity is high. We detect the narrow ra- 
diative recombination continua of O viii and O vii at their 
respective absorption edges. This means that the line ex- 
citation mechanism in EXO 0748—67 is via radiative cas- 
cades following photoionization. 

While the intensity of the continuum level varies sig- 
nificantly, the shape of the continuum is constant. In addi- 
tion, the equivalent widths of the spectral lines vary with 
the continuum intensity, but the actual emission line fluxes 
remain constant thoughout the observations. The absorp- 
tion structure is most prominent when the continuum in- 
tensity is high, but the optical depths at the edges do not 
change. We must therefore always be looking through the 
same material. 



3.3. Velocity Broadening 

All of the emission lines show significant velocity broad- 
ening. This is illustrated in Fig. ^, which shows the O Viii 
Lya line overlaid with the instrument line spread func- 
tion (LSF). Fitting the Lya lines by convolving the LSF 
with gaussian distributions, we measure velocity widths of 
CT = (2600 ±490) km s-i for Nex Lya, (1390 ± 80) km s^^ 
for Oviii Lya, and (850 ± 180) km s^^ for Nvii Lya. We 
find a direct correlation between the magnitude of the ve- 
locity width and the degree of ionization. We measure the 
shifts in the line centroids to be less than 0.020 A for all 
of the lines. This means that the systemic velocity of the 
line emitting plasma with respect to the line of sight is 
less than 300kms~^. 
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Fig. 2. Spectra for the periods of fow emission, active vari- 
ation and burst as defined in Sect. 3.1. The cumulative 
exposure time for each spectrum is indicated. 
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Fig. 3. O VIII Lya. The histogrammed data is overlaid 
with the normalized raytraced LSF. The error bars on 
the data contain both the statistical errors and the 10% 
uncertainty in the efl^ective area. 



3.4. Density Sensitivity 

We can use the line ratios in the He-like complexes to 
put Hmits on the density of the plasma in the different 
ionization regions. At sufficiently high densities the upper 
energy level of the forbidden line transition is coUisionally 
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Fig. 4. O VII He-like series. The instrument LSF, broad- 
ened to account for a 1390 km s~^ velocity field, is over- 
laid on the histogrammed data. The contributions from 
the resonance line (r), intcrcombination lines (i), and for- 
bidden line (f) are shown with thin lines. The thick line 
shows the combined fit. 



depopulated through the upper level of the intcrcombina- 
tion lines. The ratio of the forbidden line flux to the inter- 
combination line flux is therefore sensitive to the density of 
the plasma. The He-like complexes in EXO 0748—67 show 
the weak resonance lines characteristic of a photoionized 
plasma, bright intercombination lines and no obvious for- 
bidden lines (see Fig. ^. We have measured upper-limits 
to the flux in the forbidden lines and compared them to 
the flux in the intercombination lines to determine lower- 
limits to the electron densities. For O Vii we find an upper 
limit of 0.19 for the ratio, which corresponds to a lower- 
limit to the electron density of 2.0 x 10^^ cm~'^. For Neix 
we measure an upper-limit of 0.21, which corresponds to 
a lower-limit of 7 x lO^^ ^^-3 (Porquet & Dubau ^OOOj) . 



4. Discussion 

From the observed spectral behavior and the derived 
source parameters we can construct an empirical model of 
the circumstellar material. The large equivalent widths of 
the emission lines require that the emission regions sub- 
tend a fairly large solid angle in order to be sufliciently 
illuminated by the central source. These regions must also 
extend high above the plane of the binary or else we would 
have observed eclipses or other modulations correlated 
with the orbital phase. We can estimate the geometry of 
this material by comparing the observed absorption and 
emission in the hydrogen-like and helium-like oxygen ions. 
For a plasma in photoionization equilibrium the ionization 
rate is equal to the recombination rate everywhere in the 
gas. This means that the number of photons that are ab- 
sorbed by a particular ion can be directly compared to 
the number of photons emitted by that ion. Given the 
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efficiency of cascade through a particular line transition, 
Tyiino, this can be expressed as 



Mi, 



^linc -^rccombincd — ^/linc-^absorbcd ■ 



(1) 



where N is the photon flux. The ratio of observed ab- 
sorption to observed emission is then determined by the 
geometry of the system. For example, for a disk oriented 
with its axis along the line of sight, photons will be ab- 
sorbed in the surrounding material but then re-emitted 
isotropically, so that the observed line emission would be 
higher than expected from the observed absorption. For 
EXO 0748—67 we have compared the number of photons 
absorbed at the O viii and O vii edges to the number of 
photons emitted in the O viii Lya and O vii He-like com- 
plexes. We used cascade efficiencies (44% for O viii and 
75% for O vii) calculated with the at omic model for a 
recombining plasma used in Sako et al. ( f 999| ). We see ap- 
proximately three times more absorption than expected 
in each of the ions. The absorbing material must there- 
fore be flattened along our line of sight with a solid angle 
of roughly ^47r. Since we know that the system is highly 
inclined, this material must be aligned with the plane of 
the binary and therefore with the accretion disk. 

The fact that the emission lines show large velocity 
broadening with no systemic Doppler-shifts suggests that 
the observed widths are not associated with net outflow 
or inflow but rather with orbital motion around the cen- 
tral source. Using the ionization parameters of formation 
calculated as described in Sako et al. (199£) (see Table ^ 
for these values) we find a correlation with velocity that is 
consistent with an orbital velocity structure; Nex, which 
has the largest ionization parameter and is preferentially 
emitted closest to the central source, is observed to be 
moving with the largest velocity while Nvii, which has 
the lowest ionization parameter and should be emitted 
farthest from the source, is observed to be moving with 
the lowest velocity. 

The EXO 0748—67 system must therefore contain a 
thickened accretion disk where the orbiting material ex- 
tends high above the binary plane. The plasma in the 
uppermost regions probably produces both the emission 
and absorption features in the observed spectra. The con- 
tinuum intensity is independent of orbital phase so the 
central source region must be extended; a compact X-ray 
source would show orbital modulation in such a highly in- 
clined system. Since the flux in the emission lines remains 
constant, the source luminosity must be constant; varia- 
tions in luminosity would cause changes in the ionization 
structure and therefore changes in both the line flux and 
the line ratios. The optical depths at the absorption edges 
are also constant so we must always be looking through the 
same absorbing and emitting material. The observed vari- 
ations in continuum intensity are best explained by local 
obscurations of the central source region. Interestingly, the 
velocity parameters derived for the optically-emitting ma- 
terial ( ~ 2000kms"i line broadening with 2f0±92 km s"^ 
amplitude Doppler-shifts (Crampton et al. I986|) ) are very 
similar to those measured here. Although it must be in a 



very different region of ionization, it is possible that the 
optically-emitting material is also orbiting in this thick- 
ened disk. 

As a simple test of this model we can compare the 
emission measure for an orbiting plasma to the empir- 
ical emission measure. Assuming a velocity profile with 
V = {^^y^^ and a mass of I.4Mq for the neutron star, 
we can calculate the radial distances for each ion from 
the measured velocities. Then, using the ionization pa- 
rameters of formation (^ = ■:^^^) and the radial distances 
we can calculate the densities in these regions assuming a 
source luminosity of I x lO'^^ergs s"^. Defining the emission 
measure as 



EM 



dVni 



we take Ar 



and Afi 



(2) 



2- , coxi^ — ^|47r to account for the 
estimated geometry. The resulting emission measures as 
well as the components of the calculation are given in 
Table |l|. We find emission measures that range from 
- 22 X I0"cm-3 fQj. j^gx Lya to ~ 10 x lO^^cm-^ for 
Nvii Lya. 

To calculate the empirical emission measure we divided 
the line luminosity by the line power evaluated at the ion- 
ization parameter of formation. The line power for each 
ion was calculated using the HULLAC atomic codes (Bar- 
Shalom et al. 1998| ) following the procedure described in 
Sako et al. (1999). We adjusted our measured line fluxes to 
account for interstellar absorption using the cross sections 
of Morrison & McCammon (1983) with a column density 
of 1.1 X lO^^cm"'^. We assumed a distance to the source of 
D = lOkpc (Gottwald et al. 1986). The resulting emission 



measure for each ion, assuming solar abundances (Anders 
& Grevesse 1989 ), is included in Table 0. For Nex and 
Nvii we find emission measures of 2.1 and 2.6 x lO^^cm"^. 
For Oviii we find an emission measure of 0.6 x lO'^^cm"^. 
This suggests an underabundance of oxygen that may have 
interesting implications for the evolutionary state of the 
companion star. Adopting the factor of two uncertainty in 
the empirical emission measure that is estimated by Sako 
et al. ( |l999| ), we find that the predicted emission mea- 
sures are roughly a factor of five higher than the empir- 
ical emission measure except at O Viii. Considering that 
the predicted emission measures are highly dependent on 
our geometric assumptions and on a velocity profile that 
is only strictly valid for a thin Keplerian disk, the agree- 
ment with the empirical emission measures is remarkably 
good. 

This empirical model of a thickened disk is difficult to 
understand theoretically. For a disk in hydrostatic equilib- 
rium, the vertical height at a particular radius should be 
roughly given by the product of the radius and the ratio 
of the sound speed of the gas to its local orbital velocity. 
Here, we can directly infer the sound speed from the ion- 
ization parameter and the orbital velocity from the widths 
of the lines, yet we infer a disk height that grossly violates 
this condition. However, despite the theoretical challenge 
that this presents, we find no alternative model that sat- 
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Table 1. Measured and Derived Parameters for the Lya 
emission lines 



Line 


Velocity 


Radius 




Density 


EMprcdictod 


ENIenipirical 




(kms-^) 


(cm) 




(cm-3) 


(10"^crrr^) 


(10"cm-3) 


Ne X 


2600 


2.8 X 10" 


2.1 


2.0 X 10'-' 


21.8 


2.1 


O VIII 


1390 


9.6 X 10" 


1.9 


2.6 X 10" 


15.8 


0.6 


N VII 


850 


2.6 X 10^° 


1.8 


4.8 X lO" 


10.3 


2.6 



isfies the empirical constraints. The large observed line 
equivalent widths and the clear evidence for intrinsic ab- 
sorption essentially guarantee that this material is well out 
of hydrostatic equilibrium. Hopefully, as high-resolution 
spectroscopic observations of other LMXB's accumulate, 
we may gain a clearer understanding of how such accretion 
flows can be formed and maintained. 
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